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ABSTRACT 

The  effects  of  vibrational  and  chemical  nonequilibrium  on 
turbine  performance  were  investigated  separately.   The  vibra- 
tional model  was  taken  as  a  pure  nitrogen  expansion,  and  the 
chemical  model  was  taken  as  the  combustion  products  of  a 
stoichiometric  mixture  of  kerosine  and  air.   The  loss  of  per- 
formance of  fully  frozen  flow  with  respect  to  the  equilibrium 
flow  for  each  model  was  determined.   The  extent  to  which  non- 
equilibrium  will  occur  was  investigated  within  limited  ranges 
of  pressure  and  temperature.   Vibrational  nonequilibrium  can 
result  in  losses  up  to  seven  percent  with  respect  to  equilib- 
rium flow.   Vibrational  freezing  will  be  virtually  complete  a 
short  distance  after  the  throat  of  the  first  stator  assembly. 
The  losses  due  to  chemical  nonequilibrium  are  insignificant 
compared  to  the  vibrational  nonequilibrium  losses. 
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I.   INTRODUCTION 

The  aircraft  of  tomorrow  must  be  capable  of  operating  at 
high,  Mach  number,  have  large  amounts  of  excess  thrust,  be 
capable  of  operating  at  high  altitude,  and  yet  be  very  eco- 
nomical.  As  the  requirement  for  these  aircraft  becomes  more 
prevalent,  designers  of  both  airframes  and  power  plants  must 
strive  to  match  their  respective  ideas  to  produce  a  viable 
design.   The  airframe  specialist  is  concerned  with  drag  and 
high  skin-temperature  problems  which  he  hopes  to  resolve  by 
reduction  of  frontal  area  and  by  lighter  gross  weights.   The 
power  plant  specialist  investigates  higher  turbine  inlet  tem- 
peratures because  of  the  increased  efficiency  of  the  cycle 
with  high  peak  temperatures.   The  increase  in  flight  Mach 
number  induces  a  rise  in  the  temperature  of  the  cooling  air 
available  at  the  turbine.   This  higher  temperature  reduces  the 
cooling  potential  of  the  air. 

The  next  generations  of  the  turbojet  engine  will  utilize 
nearly  double  the  present  turbine  inlet  temperature  by  burn- 
ing stoichiometrically.   As  the  turbine  inlet  temperature  is 
increased,  the  effects  of  dissociation  and  vibrational  excita- 
tion play  a  significant  role  in  the  energy  balance.   The  jump 
in  turbine  temperature  for  stoichiometric  combustion  will 
produce  almost  triple  the  degree  of  vibrational  excitation,  a 
large  increase  in  local  velocities  for  the  same  local  Mach 
number  within  the  engine,  and  a  substantial  rise  in  the  degree 
of  dissociation  of  the  products  of  combustion. 
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The  phenomenon  of  vibrational  energy  freezing  occurs  due 
to  the  relatively  large  number  of  molecular  collisions  re- 
quired to  return  the  vibrational  mode  of  internal  energy  to 
equilibrium.   Freezing  of  chemical  composition  occurs  when 
the  gas  cannot  react  fast  enough  to  maintain  equilibrium  with 
a  rapidly  changing  temperature  and  pressure.   Both  of  these 
phenomena  are  familiar  to  the  rocket  designer.   The  relatively 
low  temperatures  and  high  pressures  characteristic  to  the 
turbojet  cycle  in  the  past  have  not  caused  these  nonequilib- 
rium  effects  to  any  significant  degree. 

The  presence  of  nonequilibrium  in  the  form  of  frozen 
vibrational  energy  and/or  frozen  chemical  composition  is  very 
probable  in  the  turbine  section  of  an  advanced  high  tempera- 
ture turbojet.   The  rapid  removal  of  work  from  the  gas  in  the 
turbine  reduces  both  the  temperature  and  pressure  at  a  rapid 
rate.   The  vibrational  temperature  is  the  independent  variable 
in  the  vibrational  energy  distribution,  and  there  is  a  cor- 
responding temperature  for  both  rotation  and  translation.   At 
equilibrium,  the  rotation,  translation,  and  vibration  tempera- 
tures are  the  same.   Comparison  of  the  vibrational,  rotational, 
and  trans lational  temperatures  in  the  turbine  would  show  a  more 
rapid  decrease  in  the  latter  two  relative  to  the  rate  of 
vibrational  temperature  decrease,  primarily  due  to  the  number 
of  collisions  required  to  maintain  equilibrium.   The  magnitude 
of  the  relaxation  time  relative  to  the  short  transit  time  in 
the  turbine  suggests  that  a  significant  degree  of  vibrational 
freezing  may  result.   The  effect  of  lower  pressure  would  be  to 
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reduce  the  collision  frequency,  thereby  increasing  the  relaxa- 
tion times . 

Evaluation  of  the  loss  of  available  energy  due  to 
vibrational  freezing  becomes  rather  complicated  in  multi- 
component  gases  and  involves  complex-numerical  computer  pro- 
grams.  Stollery  and  Park  [17]  and  Wilson,  et  al.  [22]  give 
examples  of  these  programs.   Exact  calculations  of  recombina- 
tion and  dissociation  rates  of  a  mixture  of  gases  such  as 
combustion  products  also  involve  complex  computer  programs. 
The  complexity  of  the  entire  problem  is  a  result  of  the  num- 
ber of  nonlinear  equations  involved  and  the  lack  of  an  exact 
solution  to  the  flow  problem  in  rotating  machines. 

This  thesis  attempts  to  apply  the  nonequilibrium  analysis 
of  vibrational  and  chemical  freezing  to  the  turbine  section 
of  a  turbojet  cycle.   The  presence  and  degree  of  each  of  these 
problems  are  investigated,  and  the  effects  of  changing  the 
thermodynamic  parameters  of  the  cycle  are  discussed. 
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II.   ANALYSIS 

High,  velocity  flows  frequently  create  rapid  changes  in 
the  flow  variables:   temperature,  pressure,  and  density.   The 
internal  modes  of  freedom  of  the  gas  molecules  do  not  always 
remain  in  equilibrium  or  return  to  equilibrium  fast  enough  to 
keep  pace  with  the  rapidly  changing  macroscopic  environment. 
The  return  to  equilibrium  of  these  internal  modes  is  governed 
by  internal  kinetic  equations.   The  existence  of  chemical  re- 
actions in  the  flow  introduces  the  chemical  kinetic  equations. 
There  are  characteristic  times  associated  with  each  of  these 
processes:   relaxation  time  with  the  internal  kinetic  equa- 
tions and  reaction  times  with  the  chemical  kinetic  equations. 

Maintenance  of  equilibrium  of  chemical  reaction  and  of 
internal  degrees  of  freedom  requires  a  number  of  molecular 
collisions.   It  is  obvious  that  true  equilibrium  flow  cannot 
exist.   The  average  time  for  a  unit  mass  of  gas  to  transit  a 
turbine  is  a  prime  factor  in  determining  whether  or  not 
freezing  occurs.   The  conversion  of  available  energy  to  useful 
work  by  the  turbine  is  accomplished  in  a  small  distance,  and 
present  turbine  technology  involves  Mach  numbers  in  the 
transonic  range.   These  facts,  combined  with  high  local  tem- 
peratures, produce  a  local  fluid  velocity  relative  to  the 
blade  of  the  order  of  two  to  three  thousand  feet  per  second. 
The  short  physical  length  of  the  turbine  results  in  an  ex- 
tremely low  transit  time.   If  the  characteristic  times  of 
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progression  of  the  chemical  and  molecular  phenomena  are  very 
small  compared  to  the  time  intervals  of  interest  in  the  flow, 
the  flow  may  be  considered  to  be  in  local  equilibrium.   Con- 
versely, if  these  characteristic  times  are  much  larger  than 
the  time  intervals  of  interest  in  the  flow,  the  phenomena 
may  be  considered  frozen. 

Theoretical  description  of  a  nonequilibrium  flow  requires 
incorporation  of  the  kinetic  equations  describing  the  processes 
into  the  usual  gas-dynamic  conservation  equations.   Analysis 
and  solution  of  such  systems  of  equations  are  difficult;  in 
all  but  the  simplest  problems,  numerical  methods  are  required. 
Bray  and  Appleton  [5]  and  Anderson,  et  al.  [2],  have  presented 
detailed  analytical  solutions  to  these  problems. 

Extensive  research  into  nonequilibrium  flows  was  inspired 
by  the  advent  of  space  exploration.   The  flow  around  re-entry 
vehicles  at  high  altitudes  involves  complex  nonequilibrium 
chemistry  [6].   In  rocket  propulsion  systems,  the  problems  of 
large  expansion  ratios  and  low  exit  pressures  tended  to  magnify 
the  situation.   Excellent  literature  surveys  by  Read  [15],  Rich 
and  Treanor  [16],  and  Taylor  and  Bitterman  [18]  provide  a 
myriad  of  data  on  vibrational  relaxation.   Gas  dynamic  lasers 
employ  intentional  vibrational  nonequilibrium  to  obtain  popula- 
tion inversions  [3] . 

The  complexity  of  solution  of  the  expansion  of  combustion 
products  through  a  turbine  should  be  readily  apparent.   The 
problem  must  therefore  be  simplified  in  order  to  obtain 
meaningful  conclusions  as  to  the  effects  of  nonequilibrium 
phenomena  on  the  expansion.   The  interaction  of  vibrational 
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excitation  and  chemical  reaction  is  not  solved  here,  but  some 
insight  into  the  level  of  existence  of  each  phenomenon  is 
presented. 

The  effects  and  mechanics  of  vibrational  nonequilibrium 
flow  are  considered  separately  from  the  effects  and  mechanics 
of  chemical  nonequilibrium  flow. 

A.   VIBRATIONAL  NONEQUILIBRIUM  FLOW 

In  order  to  investigate  the  effects  of  vibrational  non- 
equilibrium  in  a  turbine,  a  simplified  model  was  utilized. 
The  expansion  of  pure  gaseous  nitrogen  through  a  turbine  was 
considered  instead  of  the  more  complex  case  of  the  actual 
combustion  products.   Nitrogen  was  a  good  model  for  a  number 
of  reasons.   By  using  a  pure  gas,  no  chemical  reaction,  other 
than  dissociation,  was  present.   The  degree  of  dissociation 
in  the  temperature  range  considered  was  negligible.   The  prob- 
lem of  vibrational  freezing  could  therefore  be  isolated,  and 
the  calculation  of  thermodynamic  data  for  the  flow  was  greatly 
simplified.   In  addition,  the  mole  fraction  of  nitrogen  is 
large  in  the  gas  mixture. 

The  following  assumptions  were  expeditious: 

1.  Complete  (microscopic  and  macroscopic)  equilib- 
rium existed  at  the  exit  of  the  combustion 
chamber. 

2.  The  gas  behaved  as  an  ideal,  inviscid  fluid. 

3.  The  turbine  expansion  was  adiabatic. 

4.  The  flow  was  considered  as  quasi  one-dimensional. 

5.  The  velocity  of  the  gas  at  the  combustion  chamber 
outlet  was  considered  negligible  compared  to  the 
gas  velocity  after  the  turbine  exit. 
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A  turbine  schematic  (Figure  1)  was  defined.   At  station 
one,  corresponding  to  burner  outlet,  the  nitrogen  was  con- 
sidered to  be  in  complete  equilibrium.   Between  stations  one 
and  two,  the  gas  was  expanded  isentropically  through  the 
turbine.   The  work  extracted  by  the  turbine  was  assumed  to  be 
drawn  strictly  from  the  trans lational  and  rotational  modes  of 
internal  energy.   These  modes  were  considered  to  be  in  equi- 
librium with  each  other  during  the  entire  process  due  to  the 
relatively  small  number  of  collisions  required  for  equilibrium. 
This  model  assumes  the  independence  of  the  separate  modes  of 
internal  energy.   The  validity  of  this  procedure  may  not  be 
universal  in  its  application  but  was  utilized  here  for  sim- 
plicity of  analysis.   The  vibrational  energy  was  considered 
frozen  in  this  expansion  process.   (See  Figure  2.) 

Between  station  two  and  station  three,  the  gas  was  allowed 
to  return  to  complete  equilibrium  in  a  constant  area  duct. 
This  was  accomplished  by  the  redistribution  of  excess  vibra- 
tional energy  to  the  rotational-transational  gas  through  a 
process  equivalent  to  Rayleigh  heating.   The  rotational  and 
trans lational  modes  were  considered  to  be  always  in  complete 
equilibrium. 

The  extraction  of  work  results  in  a  decrease  of  total 

enthalpy.   A  hypothetical  stagnation  temperature  of  the 

rotational-trans lational  gas  at  station  two,  T'  ,  was  chosen 

r2 

to  represent  a  given  amount  of  work  extracted  by  the  turbine. 
The  prime  represents  the  state  of  the  rotational-translational 
gas.   The  ratio  of  specific  heats,  y,    was  taken  as  a  constant 
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and  equal  to  that  obtained  by  classical  equipartition  con- 
cepts for  a  diatomic  gas  with  only  rotational  and  transla- 
tional  modes.   This  value  is  y   =  1.4.   The  selection  of  a 
burner  outlet  stagnation  pressure  and  temperature  and  a  Mach 
number  at  station  two,  based  on  the  chosen  rotational - 
translational  CRT)  gas  temperature,  enabled  the  determination 
of  the  station  two  stagnation  pressure  through  the  use  of  the 
ideal  gas  isentropic  relationship  in  equation  (1) . 


P  ■    =   P 
t       t 
^2        1 


ml    "1  Y 


Tt 


Y-l 

(1) 


Use  of  an  isentropic  procedure  is  consistent  with  the  desire 
to  isolate  vibrational  nonequilibrium. 

Downstream  of  station  two,  the  gas  was  subjected  to  a 
Rayleigh  heating  process  due  to  the  relaxation  of  the  vibra- 
tional mode. 

2  2 

U2  U3 

c>   +  Cp    TI  +  =£  =   Cp    T'  +  ^  (2) 

RT        g        RT        g 
Utilizing  the  definition  of  stagnation  enthalpy, 

Qr  +  H-   =  H-  (3) 

where  Q   represents  the  energy  added  to  the  RT  gas  from  the 

vibrational  mode. 

Equation  (3)  may  be  rewritten  as: 

Q    =   H'  -  Hi   =   Cp   (T;   -  T;  )  (4) 

The  ideal  gas  isentropic  relationship  was  used  to  find  static 
temperature,  T',  which  was  needed  in  the  iteration  procedure. 
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Tt 

2 

T2   =   Y^I 2  (5) 

1+VM2 

The  determination  of  the  Mach  number  at  point  three  was  then 

required.   This  value  was  determined  by  an  iterative  process 

matching  the  Rayleigh  heat  added  to  the  classic  definition  of 

enthalpy.   The  enthalpy  was  considered  as 

r 

h(T)   =      Cp  dT  (6) 

^T 
o 

Since  the  internal  energy  modes  were  considered  separable,  the 

specific  heat  may  be  expressed  as 

C    =   C    +  C    +  C     =   C     +  C  (7) 

R      T      v        RT      v 


Using  equation  (7) ,  equation  (6)  may  be  expressed  as 

dT  (8) 


h(T)   =   \   C     dT  +   ]  C 

Jm         RT       J„         v 


T  T 

o  o 

The  energy  added  to  the  equilibrium  RT  gas  in  the  Rayleigh 
heating  model  must  be  identically  equal  to  the  difference  in 
the  vibrational  energy  of  the  complete  gas  between  stations 
one  and  three.   This  difference  was  taken  thusly 


rTi  rTi 

hCT^)       =       j  Cp       dT    +        \ 


T  -T  RT 

o  o 


Cp        dT  (9) 


Therefore 


11  =  r i{ 


MTn)  Cp       dT    +    Cp       CT1-TQ) 

m  V  RT 

o 
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Similarly 


h(T    j 


T 

■5  * 


dT   +   C„       (T--T    ) 

PRT      3      ° 


(10) 


Taking   the    difference   between   equations    (9)    and    (10) 


hCTx)    -  h(T3) 


T1 

=        ]        C         dT    +    C         (T,-T 
*L  v  RT 


(11) 


Now  define,  Q  ,  as  the  difference  in  the  vibrational  energy 
between  stations  two  and  three. 


Jm  v 


dT 


(12) 


Substituting  equations  (11)  into  equation  (12) 


QD   =   h(Tl)  -  h(T3)  -  ^^^3) 


(13) 


The  iteration  was  satisfied  when  q   of  equation  (4)  was  equal 
to  q  of  equation  (12) .   The  actual  solution  for  a  given  set 
of  input  parameters  was  obtained  by  entering  gas  tables  [11] 
with  the  chosen  Ml.   With  the  already  chosen  value  of  T'  ,  the 

values  of  T*  and  T*  were  obtained.  A  Mach  number  at  station 
3'  was  assumed,  and  equation  (4)  was  solved  in  the  following 
form, 


Qr   =   T*  Cp 


RT 


Tt. 


Tt. 


t 


t 


(14) 


Utilizing  information  from  gas  tables  [11]  the  temperature, 
T^f  was  also  obtained  and  equation  (12)  was  evaluated.  This 
process  was  repeated  until  Q   and  Q  were  within  0.5%  of  each 
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other.   Satisfaction  of  the  iteration  process  allowed  calcula- 
tion of  stagnation  temperature  and  pressure  at  station  three. 

The  investigation  of  performance  lost  due  to  vibrational 
freezing  required  a  reference  quantity.   This  reference  quan- 
tity was  taken  as  the  turbine  work  produced  by  the  expansion 
of  the  gas  with  full  equilibrium  between  the  modes  of  vibra- 
tion, rotation,  and  translation.   For  both  equilibrium  and 
nonequilibrium  cases,  the  turbine  pressure  ratios  were  assumed 
identical. 

Efficiency  was  defined  as  the  ratio  of  actual  work  to 

ideal  work  where  the  actual  work  was  considered  as  the  work 

obtained  from  the  vibrationally  frozen  model  and  the  ideal 

work  as  that  obtained  from  an  equilibrium  gas  expanding  to  the 

same  stagnation  pressure. 

vibrationally  frozen  work  , , c x 

n   =   >  n  .. — . * ; (15) 

v     equilibrium  work 

H(T   )  -  H(T   ) 
^1        ^3 

n 


V      H(T   )  -  H(T*  )  (16) 

Cl        C3 

where  T*   is  the  ideal  stagnation  temperature  at  the  stagnation 

pressure  of  station  three. 

The  stagnation  temperature  for  an  ideal  turbine  evaluated 
at  the  stagnation  pressure  of  station  three  was  determined  by 
isentropic  relationships.   By  definition, 

S(T.  ,  P.  )   =   S(T*  ,  P*  )  (17) 

^1    Cl  c3    C3 

and 

S(T,  p)   =   S°CT)  -  R  Jin  (P)  (18) 
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where  S° (T)  is  the  entropy  of  the  gas  at  a  given  temperature 
and  a  pressure  of  one  atmosphere.   Using  equation  (18)  in 
equation  (17) 

S° (T   )  -  R  In   P     =   S° (T*  )  -  R  In    P*        (19) 
*1  rl  ^3  t3 

It  was  chosen  that  P*   =  P*   and  the  values  of  T^   and  P,_ 

t3  t3  fcl       tl 

were  assumed.   Therefore,  an  equation  for  the  entropy  at  the 
desired  temperature  was  obtained. 


S° (T*  )   =   S°(T.  )  +  R  An 
t3  fcl 


p* 


\ 


(20) 


Solution  of  this  equation  yielded  T*  which  was  utilized  to 

t3 

compute  values  of  efficiency  from  equation  (16) .   This  effi- 
ciency represents  a  loss  of  performance  resulting  from  non- 
equilibrium  flow. 

A  method  to  identify  the  significance  of  vibrational  non- 
equilibrium  had  been  established,  and  the  question  of  actual 
presence  of  this  phenomenon  in  a  high  temperature  turbine 
remained. 

The  complexity  of  the  problem  of  vibrational  relaxation 
precludes  actual  solution  of  the  theoretical  equations  within 
the  scope  of  this  thesis,  and  experimental  results  of  other 
researchers  will  be  utilized  in  the  ultimate  discussion  of 
the  effects  of  this  phenomenon  on  turbine  flow.   For  complete 
understanding  it  is  prudent,  however,  to  develop  the  princi- 
ples involved. 

The  nitrogen  gas  was  considered  to  be  a  system  of  harmonic 
oscillators.  This  system  of  diatomic  molecules  was  assumed  to 
be  capable  of  collisional  interchange  of  energy  with  the 
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translational  and  rotational  degrees  of  freedom  of  a  heat 
bath  at  temperature  T.   The  heat  bath  was  provided  by  the 
translational  and  rotational  degrees  of  freedom  of  the  pure 
gas.   At  equilibrium,  the  distribution  of  oscillators  over 
the  energy  states  was  assumed  to  be  given  by  the  Boltzmann 
distribution. 

The  majority  of  data  on  vibrational  nonequilibrium,  both 
theoretical  and  experimental,  has  been  developed  for  the  lag 
of  vibrational  equilibrium  when  vibrational  degrees  of  free- 
dom are  being  excited.   Rich  and  Treanor  [16]  point  out  that 
there  is  evidence  that  the  relaxation  time  for  de-exitation 
may  be  one-fifth  or  smaller  relative  to  the  exitation  relaxa- 
tion time.   The  relaxation  time  is  a  measure  of  the  rate  of 
return  to  equilibrium  after  a  disturbance. 

The  equation  for  the  time  rate  of  change  of  the  number  of 
molecules  in  a  particular  energy  state  may  be  developed  using 
statistical  thermodynamics  and  quantum  mechanics.   The  rate 
of  transition  from  one  state  to  another  is  dependent  upon  the 
number  of  collisions  per  unit  time  between  the  oscillators 
and  the  molecules  of  the  heat  bath  that  possess  sufficient 
energy  to  cause  transition  and  on  the  probability  that  a 
single  collision  will  result  in  transition. 

For  gas  dynamic  purposes  the  total  vibrational  energy  is 
more  useful  than  the  molecular  energy  level  population  in 
each  energy  level.   This  total  energy  is  a  summation  over  all 
possible  levels  of  the  product  of  the  energy  in  a  particular 
level  and  the  number  of  molecules  in  that  level.   Manipulation 
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of  the  rate  of  change  of  the  molecular  level  population  equa- 
tion produces  an  equation  for  the  rate  of  change  of  total 
vibrational  energy.   Vincenti  and  Kruger  [20]  express  this  as 

dE        E*  -  E 

dt  x  UJJ 

where  E   is  total  vibrational  energy,  E*  is  total  vibrational 

energy  at  equilibrium,  and  x  is  the  relaxation  time. 

This  equation  shows  that  the  vibrational  energy  of  the 
system  will  always  move  toward  equilibrium  and  that  the  rate 
of  this  movement  at  any  instant  is  linearly  related  to  the 
degree  of  departure  from  equilibrium  at  that  instant.   Even 
though  this  equation  assumed  a  Boltzmann  distribution  only  at 
equilibrium  and  no  restriction  is  placed  on  the  amount  of 
departure  of  vibrational  energy  from  equilibrium,  Vincenti  and 
Kruger  [20]  limit  the  equation  to  small  departures  from  equi- 
librium.  The  limitation  to  small  departures  is  a  result  of 
the  anharmonicity  of  actual  molecules. 

The  relaxation  time,  x f    is  a  function  of  the  temperature 
and  pressure  of  the  heat  bath  and  of  the  characteristics  of 
the  vibrating  species.   For  practical  purposes  it  was  assumed 
that  equation  (21)  was  valid  for  any  number  of  excited  mole- 
cules and  could  be  applied  to  a  time  varying  heat  bath  as  long 

as  E*  and  x  were  taken  as  functions  of  the  instantaneous  state 
v 

of  the  heat  bath.   Under  these  assumptions,  x  should  be  tech- 
nically referenced  as  a  local  relaxation  time  where  numerical 
integration  could  produce  a  solution  for  a  particular  time 
span. 
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For  sufficiently  low  temperatures  the  original  Landau  and 
Teller  Theory  [20]  may  be  approximated  by: 

K   1/3 
x   =   |exp(^-)  (22) 

where  K  is  a  constant  dependent  upon  the  physical  properties 
of  the  species,  and  C  is  an  empirical  constant  dependent  upon 
the  type  of  heat  bath  molecules. 

In  the  present  theoretical  state,  the  evaluation  of  x  must 
be  accomplished  by  experimental  measurements  interpreted  in 
the  light  of  the  foregoing  theoretical  results.   The  methods 
of  ultrasonic  absorption,  shock  tube  inter ferome try ,  spectro- 
scopic data,  and  pressure  measurement,  have  been  used  experi- 
mentally to  produce  the  values  of  the  constants  for  use  in 
equation  C22) .   These  values  are  presented  in  Table  II. 

Taylor  and  Bitterman  [18]  present  numerous  graphs  of  both 
T  and  the  probability  per  collision  of  an  exchange  of  energy 
between  the  trans lational  and  vibrational  modes  of  a  molecule. 

Anderson,  et  al.  [2] ,  presents  correlations  between  re- 
laxation times,  temperature  and  pressure  for  nitrogen  and 
various  collision  partners.   These  correlations  were  utilized 
extensively  and  are  therefore  presented  here. 

log10(TP)N  _N   =   93  (T"1/3)  -  4.61  (23) 

1O^10(TP)N2-C02      =      10^10(TP)N2-N2  (24) 

log10(TP)N    _H   0      =      15.4    (T"1/3)    +    0.722  (25) 
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Bray  [4]  suggests  that  once  an  appreciable  deviation  from 
equilibrium  has  occurred  that  return  to  equilibrium  in  a 
nozzle  is  unlikely.   Therefore,  the  investigation  was  cen- 
tered on  the  stator  nozzle  block  of  the  turbine.   A  stator 
passage  with  area  change  from  throat  to  exit  given  by  equation 
(26)  and  shown  in  Figure  3  was  chosen, 

A/A*   =   0.15x  -  0.0 lx2  (26) 

where  A*  was  the  area  at  the  throat  and  x  was  the  distance 
along  the  mid-span  streamline  measured  from  the  throat.   While 
the  choice  of  this  particular  formula  for  change  of  passage 
area  was  arbitrary  and  not  particularly  critical,  this  choice 
is  representative  of  modern  turbines. 

Knowledge  of  the  area  change  and  selection  of  a  represen- 
tative y  enabled  determination  of  Mach  number,  and  thermo- 
dynamic data  from  gas  tables  as  presented  in  Table  III. 
Selection  of  a  total  pressure  and  temperature  at  the  burner 
outlet  permitted  evaluation  of  vibrational  relaxation  times 
at  the  throat  of  the  stator  with  area  change  given  by  equation 
(26)  for  nitrogen  to  nitrogen  collisions  using  equations  (22) , 
(2  3)  ,  or  the  graphs  of  Taylor  and  Bitterman  [18]  .   (See  Table 
IV.) 

As  stated  previously,  the  calculated  relaxation  times  had 
to  be  compared  to  a  gas  transit  time.   The  gas  transit  (dwell) 
times   were  determined  by  dividing  a  characteristic  length  by 
the  appropriate  velocity.   Under  this  concept,  a  total  transit 
time  was  obtained  by  dividing  the  distance  from  point  A  to 
point  B  along  the  mid-span  streamline  by  the  mean  velocity 
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between  those  two  points.   A  local  dwell  time  at  a  data  point 
was  obtained  by  dividing  the  distance  between  data  points  by 
the  velocity  at  the  upstream  point. 

Numerical  evaluation  of  equation  (21)  was  performed 
utilizing  calculated  relaxation  times  and  the  method  of  Hurle, 
Russo,  and  Hall  [9].   The  method  is  illustrated  in  Appendix  C. 
This  evaluation  illustrated  the  actual  amount  of  vibrational 
energy  as  compared  to  the  equilibrium  amount. 

An  aircraft  gas  turbine  using  hydrocarbon  fuels  will  have 
many  molecular  species  in  the  combustion  products.   It  is 
known  that  some  of  these  species  are  more  effective  for  vibra- 
tional de -excitation.   Consequently,  a  composite  relaxation 
time  involving  several  molecules  was  calculated.   The  molecular 
species  in  the  combustion  products  were  screened  to  ascertain 
the  molecules  most  effective  for  vibrational  de-excitation. 
This  screening  indicated  H„0  and  C0~  should  be  included.   A 
relaxation  time  was  thus  calculated  for  nitrogen  with  three 
different  collision  partners,  N„,  C0~  ,  H_0,  using  equations 
(23) ,  (24) ,  and  (25)  respectively  (Table  V) .   The  total  rate 
of  change  of  vibrational  energy  of  the  nitrogen  was  considered 
to  be  the  sum  of  three  distinct  quantities  corresponding  to 
the  change  induced  by  each  of  the  collision  partners. 

dE  dE  dE  dE 

XT' =  A^ +   B  XT' +   C^ (27) 

attotal        ^2"N2        ^2"C02        ^2~H2° 

Each  of  the  terms  on  the  right-hand  side  of  this  equation  was 
expressed  in  the  form  of  equation  (21)  and  substituted  to  form 
equation  (.2  8) 
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dE  E*-E        E*-E         E*-E 

V  V   V  V   V  V   V  , 

—I =   a  -^ —  +  B  — +  C  — —  (28) 

total        N2~N2      N2~C02      Tn2~H2° 

The  coefficients  (A,  B,  C)  were  taken  to  be  the  mole  fractions 
of  the  collision  partner.   Since  the  relaxation  times  with  N2 
and  CO-  are  equal,  equation  (2  8)  may  be  expressed  as 

dE. 


=   (E*-E  ) 


dt        v   v' 


1  "  XH20    XH20 


TN2-N2     TN2-H20 


(29) 


where  X^  0  is  the  mole  fraction  of  water  in  a  gas  mixture. 

Comparing  equation  (29)  with  equation  (21) 

,  1   "  *H  0   XH_0 

+  — - (30a) 


X  .         T  T 

composite     N2~N2       N^-H^O 

The  extension  of  equation  (30a)  to  systems  with  a  larger  num- 
ber of  components  is  obvious.   Note  that  E   is  the  vibrational 
energy  of  nitrogen  only.   Values  of  x  composite  are  in  Table 
VI. 

Consider  a  large  number  of  N*,  where  N*  denotes  vibra- 
tionally  excited  N~.   If  the  number  of  N*  exceeds  the  equilib- 
rium value,  there  will  be  a  decrease  in  the  number  of  N*  as 
time  progresses.   The  N*  may  be  de-excited  by  collisions  with 
H?0  or  COp  or  N^ .   The  collision  frequency  with  H.O,  C02 ,  or 
N?  depends  upon  the  mole  fraction.   The  probability  that  a 
given  collision  will  result  in  de-excitation  depends  on  the 
partners;  i.e.  on  C  and  K_  in  equation  (22).   Using  this 
physical  picture,  equations  (27) ,  (2  8) ,  and  (29)  follow.   In 
equation  (29) ,  each  x  is  evaluated  at  a  pressure  equal  to  that 
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exerted  by  the  gas  mixture.   Alternately,  an  expression  for 
x  composite  would  be 

1  i        1 

+  -^ (30b) 


composite      N?~N2    TNp-H-0 
where  partial  pressures  are  used  to  evaluate  t  for  each 
species.   The  analogy  of  £he  equation  for  t  composite  to  the 
formula  for  resistances  in  parallel  is  apparent. 

B.   CHEMICAL  NONEQUILIBRIUM  FLOW 

In  order  to  investigate  the  effects  and  mechanics  of 
chemical  nonequilibrium  on  the  flow  in  a  turbine,  the  products 
of  combustion  of  a  stoichiometric  reaction  between  air  and  a 
kerosine-type  hydrocarbon  fuel  were  considered  and  a  turbine 
expansion  calculation  performed.   The  choice  of  a  stoichio- 
metric mixture  ratio  was  indicated  by  the  current  trends  and 
research  goals  in  the  turbojet  industry.   The  internal  degrees 
of  freedom  of  all  components  of  the  gas  were  assumed  to  be  in 
equilibrium  throughout  the  expansion.   Even  though  there  is  a 
known  relationship  between  vibrational  modes  and  chemical  re- 
action, this  assumption  was  necessary  for  analytical  purposes. 

The  gas  was  taken  as  an  inviscid,  dif fusionless  fluid  with 
quasi  one -dimensional  expansion  through  an  ideal  turbine. 
The  process  was  considered  to  be  adiabatic.   The  gas  was  ex- 
panded first  as  chemically  frozen  and  then  in  full  chemical 
equilibrium.   The  stations  in  the  process  were  as  illustrated 
in  Figure  8.   The  concepts  used  in  this  section  closely  follow 
Penner's  [13]  development  for  chemically  frozen  flow.. 
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In   chemically    frozen    flow,    the    average   molecular  weight 
obviously   remains    constant;    therefore,    the   work    derived   from 
this   expansion    can  be   expressed   as 


Workfrozen      =      J  *i  (Hi  <%>    '  Hi  <%>  > 


(31) 


where  X.  is  the  molar  fraction  of  species  i  at  station  1,  n  is 

the  number  of  species  present,  and  H. (T)  is  the  total  molar 

enthalpy  of  species  i  at  temperature  T. 

Given  values  of  P.  ,    T.  ,  and  T   ,  the  total  pressure  at 

Cl    Cl       Z2 

station  two  was  determined  from  the  isentropic  relationship 


=   P, 


Y' 


Y*-l 


(32) 


where   y*    is    an    appropriately    defined   average    specific  heat 
ratio   during  expansion.      Penner    [13]    gives    this    ratio   as: 


n 


*      = 


I    X    (S°(T.     )     -    S°(T.     )) 
i=l    x  ^1  c2 


n 


I    Xi(S°(Tt    )     -    S°(Tt    ))     -    R   £n 


Tt 


(33) 


i=l 


The  chemical  nonequilibrium  efficiency,  n  r    was  defined 
as  the  ratio  of  the  work  extracted  from  chemically  frozen  gas 
to  the  work  extracted  from  equilibrium  gas  expanded  to  the 
same  total  pressure.   This  efficiency  is 


n, 


Hl  -H2 

Mi 

Hl  "  H3  ~ 
M 

X3 


(34) 
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where 


n 

H        =         I    X    (T       ,    P      )H    (T      )  ■       (35) 

i  i=1   i      t^         x:x      J-      t1 


n 
R        =         I    X.  (T       f    P       )H.  (T.     )  (36) 

Z  i=l    r      Cl         Cl      X      C2 


H3      =      I/i'V    V^V  (37) 

i=l  3  3  3 


n 
M,  __=        I    X.  (T       ,    P      )m.  (38) 

i  i=l   ^      rl         rl      1 


M3      -       JX.(T  P      )m.  (39) 

1=1  3  3 

Once  the  quantities  in  the  numerator  of  equation  (34) 
were  established  as  functions  of  known  or  given  quantities , 
attention  was  focused  on  the  denominator. 

For  equilibrium  flow,  the  analysis  was  restricted  to  a 
given  weight  of  gas  mixture  due  to  the  temperature  dependence 
of  the  average  molecular  weight  of  the  mixture.   For  purposes 
of  comparison  it  was  convenient  to  choose  the  initial  mole- 
cular weight  of  the  gas  mixtures  as  the  characteristic 
weight;  thus  the  ideal  or  equilibrium  flow  work  was: 

Ml 
Work    .   =   H, ±-   H,  (40) 

eqUl         1     rr     3 

m3 
where 

n  r  1 

H    =    I    X. (T.  ,  Ptn)  AH°    +  H. (T.  )  -  H. (T)      (41) 
1      ±=1  i   tx     1  |   f.      i  t±  i   oj 
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and 

«3  -  ?  xi(Tt  -  pt '  Kf. +  Hi(Tt  i  -  Hi<To>]  <42> 

1  =  1      3     3  L_    i         1  J 

At  this  point,  all  the  necessary  quantities  for  evaluation  of 

r\      were  known  except  T   which  could  be  determined  from 
C  fc3 

isentropic  considerations. 

Isentropic  flow  occurred  for  the  total  weight  of  the  gas 
mixture.   Since  M_  was  chosen  as  the  fixed  weight  for  analysis, 
the  isentropic  relationship  was: 

Rl  - 
S(T   ,  P   )   =   -i-  S(T   ,  P   )  (43) 

^1     1       M3     ^3     3 

S    (T,P)    is    the   molar  entropy   of   an   ideal    gas   mixture    at    tem- 
perature  T    and  pressure   P    and   may  be   expressed   as 

n  n 

S(T,P)       ="       I    X.(T,P)S?(T)     -    RUn    P   +      J    X. (T,P) £nX. (T,P) ) 

i=l    L  1  i=l   L  1 

(44) 

where  S°  (T)  is  the  molar  entropy  at  temperature  T  and  one 

atmosphere  of  pressure.   Equations  (4  4)  and  (44)  were  combined 

to  yield 

I    X.  (T   ,  P   )fs?(T   )  -  R  In   P    -  Rta  X.  (T   ,  P   )1 
i=l      1     1  L     1  1  1     1  J 


M,   n  r 

=      -±      I    X    (T   ,P   )  S?  (T   )  -  R  In   P 
M3  i=l  X      C3   l3  L      3  . 


-  R  In   X.  (T   ,  P   ) 

1  Z3        t3 


(45) 


This  equation  was  used  to  solve  for  T    by  iteration.   Once 

t3 

this  has  been  accomplished,  evaluation  of  equation  (34)  was 
possible. 
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The  mechanics  of  the  chemically  reacting  flow  were  then 
investigated.   In  the  combustion  chamber  of  a  turbojet,  a 
significant  portion  of  the  chemical  energy  released  is  used 
to  break,  chemical  bonds  producing  more  elementary  species. 
This  energy  is  unavailable  to  do  work  unless  it  is  first  re- 
leased through  a  recombination  of  the  dissociated  combustion 
products.   The  rate  of  each  chemical  reaction  must  be  compared 
to  a  dwell  time  of  the  flow  at  any  particular  stage  in  order 
to  determine  the  extent  of  the  unavailability  of  energy. 

Consider  the  typical  three  body  reaction 

rf 
A+B  +  M+AB  +  M  (46) 

rb 
where  M  is  a  third  body  which  removes  energy  during  the  re- 
combination and  r,.  and  r,  are  reaction  rates  in  the  forward 

f      b 

(recombination)  and  backward  (dissociation)  directions 
respectively.   This  equation  represents  two  separate  reactions 
with  the  connection  that  rf  must  equal  r,  at  equilibrium.   As 
the  temperature  and  pressure  move  toward  lower  values  during 
the  expansion  process,  the  reaction  in  equation  (46)  should 
move  in  the  forward  direction.   If  the  reaction  rate  is  not 
fast  enough  to  maintain  equilibrium  conditions,  recombination 
will  lag,  and  the  gas  driving  the  turbine  will  have  dissociated 
products  whose  dissociation  energy  has  not  contributed  to  the 
kinetic  energy  of  the  flow.   The  lag  of  recombination  results 
in  a  reduction  of  performance. 

Equilibrium  flow  calculations  involve  the  simultaneous 
solution  of  the  conservation  equations  for  mass,  energy, 
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momentum,  and  species,  plus  the  equation  of  state.   This  sys- 
tem is  simplified  for  frozen  flow  because  the  species  composi- 
tion is  constant.   Computer  programs  for  solution  of 
equilibrium  and  nonequilibrium  flows  are  presented  by  Gordon 
and  Zeleznik  [8]  and  Zupnik,  et  al.r  [23]. 

For  relaxing  flow  in  a  turbine  the  problem  is  complex 
because  the  nonlinear  kinetic  differential  equations  of 
species  concentration  and  the  conservation  of  atom  equations 
are  now  involved.   Solution  of  the  chemical  relaxation  problem 
is  obviously  desirable  but  beyond  the  scope  of  the  present  in- 
vestigation.  Penner  [12]  and  Bray  [4]  have  proposed  criteria 
to  determine  when  the  shift  may  be  made  from  equilibrium  to 
nonequilibrium  analysis.   Application  of  these  criteria  by 
various  researchers  has  produced  varied  degrees  of  success. 

Newhall  [10]  presents  twenty -eight  elementary  reactions 
(Table  X)  that  may  be  considered  significant  in  the  flow  of 
combustion  products  resulting  from  a  mixture  of  a  hydrocarbon 
fuel  and  air.   These  equations  may  be  integrated  numerically 
with  the  coupled  energy  equation  to  produce  concentration- 
time  histories  for  each  species  during  expansion,  but  this 
was  beyond  the  scope  of  this  thesis. 

The  science  of  chemical  kinetics,  the  study  of  the  rates 
of  chemical  reactions,  is  complex  and  difficult,  both  experi- 
mentally and  theoretically.   A  discussion  of  both  simple  col- 
lisional  theory  and  transition  state  theory  is  given  in  Pratt 
[14]  and  emphasizes  the  heavy  reliance,  except  in  the  simplest 
reactions,  upon  empirical  data. 
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The  chemical  reaction  rate  may  be  considered  to  be  pro- 
portional to  the  total  number  of  collisions  between  reacting 
chemical  species  and  to  the  fraction  of  this  total  in  which 
sufficient  energy  is  available  for  the  reaction  to  occur.   In 
flow  systems  the  rate  of  reaction  is  also  dependent  on  the 
mass  and  energy  transport  mechanisms.   The  collision  rate  is 
dependent  upon  both  concentration  and  effective  collision 
cross  section  of  the  molecules  involved. 

The  most  widely  utilized  approach  to  the  reaction  rate 
problem  is  through  the  use  of  the  Arrhenius  equation 

K   =   A(T)  exp(-E  /RT)  (47) 

a 

where  A  is  the  pre-exponential  factor  and  E   is  the  experi- 

a 

mental  activation  energy.   The  values  of  both  A  and  E   are 
obtained  from  experimental  plots  of  In   K  vs  l/T.   Both  A  and 
E   may  be  functions  of  temperature. 

For  elementary  reactions,  such  as  those  that  make  up  the 
complex  reactions  in  the  combustion  products,  the  rate  con- 
stant, K,  from  the  Arrhenius  equation  (47)  may  be  employed  as 
follows:   The  rate  of  change  in  the  concentration  of  a  product 
AB  produced  by  a  reaction  such  as  equation  (46). 

A+B    +    M+AB    +   M  (46) 


«- 


is   equal   to 


dd^AB)       =      K    (A)     (B)     (M)  (48) 


for  an  elementary  reaction  where  the  parentheses  signify 
concentration  in  moles  per  cc. 
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Each  of  the  equations  in  Table  X  had  corresponding  for- 
ward and  backward  rate  constants  from  Table  XI.   Each  of 
these  rate  constants  was  calculated  at  the  throat  of  the 
turbine  stator  model  introduced  previously  (Table  XII) .   The 
mole  fractions  of  the  species  in  the  combustion  products  were 
calculated  at  the  stator  throat  using  Banes,  et  al.,  [3]  and 
are  presented  in  Tables  XIII  and  XIV.   A  conversion  factor 
that  enabled  the  use  of  mole  fractions  instead  of  molar  con- 
centrations in  equation  (48)  was  developed  and  is  presented  in 
Appendix  D.   The  reaction  rates  for  the  equations  considered 
were  calculated  at  the  stator  throat  and  are  tabulated  in 
Table  XV.   The  total  reaction  rate  for  each  specie  was  deter- 
mined by  a  summation  of  the  rate  of  destruction  or  creation  of 
the  specie  from  the  applicable  equation.   If  the  particular 
specie  appeared  on  the  left-hand  side  of  the  reaction  equa- 
tion in  Table  X  it  was  considered  as  being  destroyed  at  the 
forward  rate  and  produced  at  the  backward.   In  Table  XI,  the 
total  reaction  rate  at  the  stator  throat  for  each  specie  is 
presented. 

The  determination  of  the  presence  of  chemically  frozen 
flow  in  the  turbine  required  establishment  of  a  characteris- 
tic reaction  time.   The  total  reaction  rate,  R . ,  of  a  specie 

-1    -3 
i,  had  units  of  moles  sec    cm   .   Multiplying  R.  by  the 

3 

specie  molecular  weight  produced  units  of  grams  of  i  per  cm 

per  sec,  which  yields  a  rate  of  change  of  density. 

r.  m.   =   __i  (49) 

i   i      dt 

By  taking  the  gas  average  density  at  the  throat,  p".  ,  and 
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dividing  both  sides  of  equation  (49)  by  p.,  a  rate  of  change 
of  mass  fraction  is  obtained  for  small  mixture  density 
changes . 


R.  m. 

1   1. 


Pidt 


By  definition 


m. 
Y.   =   —  X. 

M 

For  sufficiently  small  changes 
dY. 


Y.       -   Y. 
X2  Xl 


dt 


fc2    "   fcl 


dY, 

_: 

dt 


(50) 


(51) 


(52) 


Letting   t„    -   t,    =    t        equation    (50)    may  be    rewritten   as 

i 


R.    m. 


dY. 

3. 

dt 


Y.       -   Y. 
fcR. 


or 


'R, 


m.  m. 

x.     =r    -    x.     -± 

12    M2  xl   M, 

r!    nT 

r      i 


(53) 


Rearranging 


R. 


R. 


X. 
l. 


-t 


M, 


M, 


(54) 


2       "1 

Evaluation  of  equation  (5  4)  provided  a  basis  for  comparing 
reaction  times  with  dwell  times. 
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III.   RESULTS  AND  DISCUSSION 

Values  of  the  vibrational  efficiency  defined  by  equation 
(16)  were  calculated  using  stagnation  conditions  of  2200 °K 
and  20  Atmospheres.   These  values  are  presented  in  Table  I. 
The  data  in  Table  I  are  plotted  in  Figure  4  as  loss  percentage 
against  the  change  of  total  temperature  of  the  RT  gas  from 
stations  1  and  2.   See  Figure  1  for  station  numbering.   Selec- 
tion of  this  abscissa  illustrates  the  variance  of  the  vibra- 
tional efficiency  with  the  amount  of  work  extracted. 

The  data  for  different  Mach  numbers  indicates  that  this 
efficiency  as  defined  is  a  definite  function  of  turbine  exit 
Mach  number.   It  is  known  that  for  a  given  static  temperature 
there  is  a  known  amount  of  energy  stored  in  the  vibrational 
mode  of  internal  energy  at  equilibrium.   As  the  exit  Mach  num- 
ber increases,  the  kinetic  energy  extracts  a  larger  portion  of 
the  rotational  and  translational  energy  depicted  in  Figure  2. 
Therefore,  the  static  temperature  of  the  RT  gas  is  lower  and 
the  amount  of  equilibrium  vibrational  energy  is  lower.   The 
total  vibrational  energy  at  the  turbine  exit  is  the  same  as 
the  equilibrium  vibrational  energy  at  the  burner  outlet.   As 
the  Mach  number  increases,  the  RT  gas  static  temperature  de- 
creases and  the  fraction  of  the  total  vibrational  energy  at 
the  exit  that  corresponds  to  the  RT  gas  temperature  decreases. 
Therefore,  the  flow  experiences  a  larger  departure  from  equi- 
librium and  an  attendant  increase  in  the  losses. 
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The  apparent  minima  of  curves  in  Figure  4  as  indicated  by 
the  dashed  line  is  a  result  of  the  mechanism  of  work  removal. 
The  work  was  extracted  from  a  RT  gas  which  had  both  modes 
fully  excited  in  the  temperature  range  considered.   This  full 
excitation  made  the  work  extraction  a  linear  function  of  tem- 
perature.  The  equilibrium  work  in  the  denominator  of  equation 
(15)  was  drawn  from  all  modes  of  the  gas.   At  the  temperatures 
considered,  vibrational  specific  heat  is  nonlinear  and  is  the 
apparent  cause  of  the  minima. 

There  were  three  sources  of  vibrational  relaxation  times: 
equation  (22) ,  the  graphs  of  Taylor  and  Bitterman  [18] ,  and 
the  equations  of  Anderson,  et  al.,  [2].   In  Table  IV,  relaxa- 
tion times  for  nitrogen-nitrogen  collisions  are  calculated  by 
the  three  methods.   The  values  are  for  an  area  change  as  de- 
fined in  equation  (26)  and  stagnation  conditions  of  2200 °K  and 
20  atmospheres.   The  variation  in  the  values  of  relaxation 
time  when  calculated  by  the  three  methods  indicates  that  the 
data  basis  for  each  calculating  method  is  different.   The 
graphs  of  Taylor  and  Bitterman  [18]  had  insufficient  range 
.for  this  analysis.   The  values  obtained  through  the  use  of 
equation  (22)  seem  unusually  small  and  therefore  were  not  used. 
The  Anderson,  et  al. ,  [2]  data  is  a  later  development  with 
information  for  many  different  collision  partners,  and  its 
use  allowed  more  flexibility  of  investigation. 

The  relaxation  times  for  nitrogen  as  a  function  of  stator 
area  ratio  are  given  in  Table  V  for  a  stagnation  temperature 
of  2200 °K  and  various  stagnation  pressure  ranging  from  5  to 
20  atmospheres.   Equation  (26)  was  used  to  calculate  area  ratio 
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as  a  function  of  distance  along  the  midspan  streamline.   These 
data  are  plotted  in  Figure  5  with  the  area  ratio  as  the 
abscissa.   The  area  ratio  values  are  of  significance  only  in 
that  they  define  data  collection  points.   The  pressure  de- 
pendence of  the  relaxation  is  aptly  demonstrated  in  this 
figure.   The  lower  pressure  results  in  a  lower  collision  fre- 
quency and  thereby  slows  the  relaxation.   The  curve  of  relaxa- 
tion time  for  nitrogen-water  collisions  emphasizes  the 
catalytic  effect  of  water  on  the  relaxation  process.   The 
presence  of  water  in  the  combustion  products  of  a  hydrocarbon 
and  air  mixture  helps  to  delay  the  onset  of  vibrational 
freezing. 

The  dwell  time  as  plotted  in  Figure  5  represents  the  re- 
ciprocal of  the  local  velocity  in  units  of  seconds  per  inch 
multiplied  by  a  factor  of  ten.   The  order  of  magnitude  of  the 
dwell  time  was  raised  because  exact  numerical  solution  of 
equation  (21)  indicated  vibrational  freezing  is  significant 
at  A/A*  =  1.14.   The  factor  of  ten  multiplication  results  in 
the  dwell  time  and  relaxation  time  curve  intersection  agree- 
ing with  that  solution.   The  multiplication  factor  makes  the 
dwell  time  longer  and  more  pessimistic.   The  intersection  of 
the  dwell  time  curve  with  the  relaxation  time  curves  represents 
a  position  downstream  of  which  the  flow  may  be  considered  to 
be  vibrationally  frozen.   For  the  pure  nitrogen,  this  point 
occurs  at  a  low  area  ratio  indicating  that  the  flow  is  prob- 
ably fully  frozen  for  practical  calculations.   Once  the  flow 
has  frozen,  Bray  [4]  indicates  that  it  will  remain  frozen  or 
at  least  in  some  form  of  nonequilibrium.   The  fact  that 
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freezing  is  indicated  early  in  the  stator  section  of  the 
turbine  implies  that  the  flow  should  be  considered  as  vibra- 
tionally  frozen  throughout  the  turbine  performance 
calculations . 

Equation  (21)  was  evaluated  for  various  conditions  and  the 
data  plotted  in  Figures  6  and  7.   Both  of  these  figures  dra- 
matically demonstrate  the  departure  from  vibrational  equilib- 
rium that  occurs  at  a  low  area  ratio.   The  effects  of  other 
collision  partners  are  shown  in  Figure  7  using  data  from 
Tables  VII  and  VIII.   The  upward  movement  of  the  curve  in 
Figure  7  when  the  midspan  streamline  length  is  reduced  is  a 
result  of  a  lower  gas  transit  time. 

Values  of  the  chemical  nonequilibrium  efficiency  were  com- 
puted from  equation  (34)  for  stagnation  conditions  of  2200°K 
and  20  atmospheres  of  pressure  at  the  burner  outlet.   These 
values  are  presented  in  Table  IX.   Data  from  Table  IX  were 
plotted  as  loss  percentage  in  Figure  9  with  the  change  of 
total  temperature  from  burner  outlet  to  turbine  exit  as  the 
abscissa.   Here  again  the  total  temperature  change  is  an  indi- 
cator of  the  amount  of  work  extracted.   The  losses  due  to 
chemical  freezing  in  this  stoichiometric-mixture  expansion  are 
rather  low.   For  the  temperature  range  considered,  the  amount 
of  nitrogen  dissociation  is  negligible.   This  is  significant 
because  nitrogen  comprises  nearly  three-quarters  of  the  gas 
mixture.   For  the  stoichiometric  mixture  there  is  very  little 
oxygen  remaining  to  dissociate  after  the  reaction.   Consider- 
ing the  size  of  the  losses,  the  tendencies  of  the  curve  are 
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subject  to  question.   The  precision  of  the  thermochemical 
data  utilized  for  the  calculation  would  have  significant  ef- 
fect on  the  results. 

The  characteristic  reaction  time  for  each  specie  was 
calculated  from  equation  (54)  using  the  data  for  the  stator 
throat  and  stator  exit  at  total  temperature  of  2200 °K  and 
total  pressures  of  15  and  5  atmospheres.   These  times  are  in 
Table  XVIII.   Only  the  reactions  considered  to  be  the  fastest 
were  utilized  in  the  analysis.   The  influence  of  slower  reac- 
tions upon  the  fast  reactions  was  not  analyzed.   Each  of  the 
times  was  analyzed  and  a  classification  was  performed.   The 
gas  transit  time  from  stator  throat  to  exit  was  compared  to 
each  specie  reaction  time,  and  the  specie  was  placed  in  an 
appropriate  section  of  Table  XIX. 

A  negative  reaction  time  for  a  specie  indicated  that  the 
predominate  or  fast  reactions  involving  that  specie  were  away 
from  equilibrium  for  the  particular  conditions  of  temperature, 
pressure  and  concentrations.   The  real  life  case,  including 
reactions  not  considered  here,  would,  for  long  time  periods, 
cause  the  reaction  times  to  become  positive.   The  negative 
reaction  time  species  were  classified  as  frozen.   A  long  or 
short  positive  reaction  time  compared  to  the  transit  time 
classified  the  specie  as  frozen  or  equilibrium  respectively. 
Specie  reaction  times  of  the  same  order  of  magnitude  as  the 
transit  times  were  classified  as  unknown.   The  approximation 
used  to  calculate  the  reaction  times  was  not  sufficiently 
accurate  to  make  a  definite  statement  as  to  the  freezing  of 
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the  particular  specie.   Table  XIX  definitely  indicates  the 
presence  of  chemically  near-frozen  flow. 

Several  aspects  of  the  current  state  of  the  art  in  turbo- 
jet development  tend  to  reinforce  the  gravity  of  the  losses 
calculated  in  this  paper.   As  the  aircraft  operate  at  higher 
and  higher  altitudes,  the  total  pressure  at  the  burner  outlet 
will  decrease.   This  will  insure  the  presence  of  vibrational 
and  chemical  nonequilibrium  as  illustrated  in  Figure  5,  6,  7, 
and  Table  XIX.   The  use  of  hydrogen  fuel  would  produce  large 
amounts  of  water  thereby  delaying  the  vibrational  freezing 
effects.   The  magnitude  of  the  burner  outlet  temperature  will 
vary  as  the  ability  of  the  turbine  blade  cooling  method  to 
handle  high  inlet  temperature,  and  this  variation  will  effect 
the  degree  of  nonequilibrium  present. 
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IV.   CONCLUSIONS 

The  vibrational  freezing  of  the  flow  in  the  turbine 
section  of  an  advanced  turbojet  using  a  stoichiometric  mix- 
ture of  kerosine  and  air  can  result  in  losses  of  up  to  seven 
percent  when  compared  to  equilibrium  flow. 

The  vibrational  freezing  of  the  flow  will  be  virtually 
complete  a  short  distance  after  the  throat  of  the  first 
stator  assembly. 

The  losses  due  to  chemical  freezing  are  insignificant 
compared  to  the  vibrational  freezing  losses. 
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APPENDIX  A 


TABLE  I 


VIBRATIONAL  EFFICIENCY  WITH  TOTAL  TEMPERATURE  OF  220 0°K 
AND  TOTAL  PRESSURE  OF  20  ATM 


T£,  °K 

M'  =  .1 

M'  =  .3 

M2  "  -5 

1500 

96.37 

95.21 

92.66 

1550 

96.85 

95.40 

92.63 

1600 

97.00 

95.52 

92.49 

1650 

97.27 

95.95 

92.43 

1700 

97.71 

95.89 

92.02 

1750 

97.90 

95.80 

91.62 

1800 

97.98 

95.80 

91.03 

1850 

97.81 

95.43 

90.51 

1900 

98.14 

95.14 

89.46 

TABLE    II 

VALUES    OF    CONSTANTS    FOR   USE    IN    EQUATION     (22) 
FROM   VINCENT  I    AND   KRUGER    [20] 


Species 

Heat  Bath 
| Molecules 

C 
|  Atm-mi  crose  c 

1  K2 

°K    1 

Approximate  Range 
of  T  °K 

°2 

°? 

5.42 

X 

io-5 

2.95 

X 

io6 

800  -  3200 

°2 

Ar 

3.5 

X 

ID'4 

2.95 

X 

io6 

1300  -  4300 

N2 

N2 

7.12 

X 

IO"3 

1.91 

X 

io6 

800  -  6000 

NO 

NO 

4.86 

X 

io"3 

1.37 

X 

io5 

1500  -  3000 

NO 

NO 

6.16 

X 

io"1 

1.37 

X 

10  5 

1500  -  4600 
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TABLE  III 

STATOR  AREA  CHANGE  WITH  ASSOCIATED  IDEAL  GAS  DATA 
FROM  WANG,  ET  AL.  [21] 
y  =  1.26 


A/A* 

M 

T/Tt 

P/Pt 

Velocity, 
M/sec 

Dwell  Time, 
sec/in 

1 

1.0 

.885 

.5531 

854 

29.7 

1.14 

1.43 

.791 

.3191 

115  4 

22 

1.27 

1.59 

.7526 

.2523 

1252 

20.3 

1.37 

1.69 

.7292 

.217 

1310 

19.4 

1.46 

1.77 

.7106 

.1925 

1354 

18.8 

1.51 

1.81 

.702 

.18 

1376 

18.5 

1.55 

1.83 

.696 

.172 

1389 

18.3 

TABLE    IV 

COMPARISON    BETWEEN    VARIOUS    METHODS    OF    OBTAINING    STATOR 
RELAXATION    TIMES    WITH    TOTAL    TEMPERATURE    =    2200 °K 
AND   TOTAL    PRESSURE    =    15    ATM 


A/A* 

T 

P 

TEqn.     i 

[22) 

Graphs 

TEqn.     (23) 

1 

1947 

11.06 

13.3 

63. 

3 

61.1 

1.14 

1740 

6.38 

33.7 

172. 

4 

207.5 

1.27 

1656 

5.05 

50.6 

238 

353 

1.37 

1604 

4.34 

65.8 

** 

498 

1.48 

1563 

3.85 

81.3 

** 

658 

1.51 

1544 

3.60 

90.8 

** 

759 

1.55 

1531 

3.44 

98 

** 

* 

837 

**    Unable    to  extrapolate    the    graph    further. 
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TABLE    V 

STATOR   RELAXATION    TIMES    IN    MICRO-SECONDS 
OF  NITROGEN    MOLECULES    UNDER   VARIED    TOTAL    PRESSURES 
WITH    TOTAL    TEMPERATURE    =    2200°K    USING   EQUATION     (23) 


a  /a  * 

Pressure,    Atm 

±\/  £\ 

20 

15 

15 

10 

5 

** 

1 

61.1 

82.9 

10.8 

12  4 

249 

1.14 

207.5 

277 

21 

415 

830 

1.27 

353 

471 

28 

706 

1412 

1.37 

498 

663 

33 

996 

1992 

1.46 

658 

876 

39 

1316 

2631 

1.51 

759 

1011 

42 

1577 

3035 

1.55 

837 

•1116 

44 

1673 

3347 

**    Collision   Partner   is   H_0   in   this    column   vice   N~ 
as    in   the    other    columns. 


TABLE    VI 

COMPOSITE    RELAXATION    TIMES    CALCULATED   USING   EQUATION     (30a) 

WITH   TOTAL    PRESSURE    =    15    ATM 


A/A* 

N2- 

■N2 

H20-N2 

Composite 

1 

82. 

9 

10.8 

44.2 

1.14 

277 

21 

106.4 

1.27 

471 

27.9 

152.5 

1.37 

663 

33.4 

190.6 

1.46 

876 

38.7 

228 

1.51 

1011 

42 

250.5 

1.55 

1116 

44.3 

267 
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TABLE    VII 

RATIO   OF   ACTUAL   VIBRATIONAL    ENERGY 
TO   EQUILIBRIUM   VIBRATIONAL    ENERGY    AT    THE    THROAT 
UNDER   VARIOUS    TOTAL    PRESSURES    AT    TOTAL    TEMPERATURE    =    2200 °K 


Total 

E* 

E 

Total 

E 

Pressure,  ATM 

V 

E* 
Vth 

A/A* 

V 

E* 
Vth 

Pressure 

A/A* 

V 

E* 
Vth 

20 

1 

1 

1 

10 

1 

1 

.8388 

1.14 

.9919 

1.14 

.9958 

.7733 

1.27 

.9779 

1.27 

.9886 

.7331 

1.37 

.9634 

1.37 

.9831 

.7006 

1.46 

.9553 

1.46 

.9786 

.6856 

1.51 

.9485 

1.51 

.9749 

.6755 

1.55 

.9423 

1.55 

.9715 

15 

1 

1 

•  5 

1 

1 

1.14 

.9938 

1.14 

.9979 

1.27 

.9  832 

1.27 

.9942 

1.37 

.9751 

1.37 

.9914 

1.46 

.9687 

1.46 

.9891 

1.51 

.9633 

1.51. 

.9871 

1.55 

.9584 

1?55 

.9854 
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TABLE  VIII 

RATIO  OF  ACTUAL  VIBRATIONAL  ENERGY  TO  EQUILIBRIUM  VITRATIONAL 
ENERGY  AT  THE  THROAT  WITH  TOTAL  PRESSURE  =  15  ATM  AND  TOTAL 
TEMPERATURE  =  2200°K  FOR  VARIOUS  MIDSPAN  STREAMLINE  LENGTHS 


E* 

V 

E* 

A/A* 

E 

V 

E* 

E 

V 

E* 

E 

V 

E* 

Vth 

Vth 

Vth 

Vth 

t 

tt 

ttt 

1 

1 

1 

1 

1 

.8388 

1.14 

.9938 

.9979 

.9946 

.7733 

1.27 

.9832 

.9942 

.9846 

.7731 

1.37 

.9751 

.9914 

.9758 

.7706 

1.46 

.9687 

.9891 

.9680 

.6856 

1.51 

.9633 

.9871 

.9609 

.6755 

1.55 

.9584 

.9854 

.9544 

t   :   Midspan  streamline  length  is  6";  collision  Partner 

V 

tt   :   Midspan  streamline  length  is  2";  collision  Partner 

N2- 

ttt   :   Midspan  streamline  length  is  2";  collision  Partners, 

N«,  C02  and  H20. 


TABLE  IX 

CHEMICAL  EFFICIENCY    CALCULATED    USING   EQUATION     (3 4] 
WITH   TOTAL   TEMPERATURE    =    2200 °K 


Turbine 
Total   Tempe 

Exit 
nature 

nc 

1500 

98.92 

1600 

99.00 

1700 

99.15 

1800 

99.15 

1900 

99.50 
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TABLE    X 
ELEMENTARY    CHEMICAL    REACTIONS    FROM   NEWHALL     [10] 


1. 

H20 

+ 

M 

OH 

+ 

H 

+ 

M 

2. 

H2 

+ 

M 

H 

+ 

H 

+ 

M 

3. 

NO 

+ 

M 

-*- 

N  . 

+ 

0 

+ 

M 

4. 

N2 

+ 

M 

->• 
-t- 

N 

+ 

N 

+ 

M 

5. 

N02 

+ 

M 

NO 

+ 

0 

+ 

M 

6. 

°2 

+ 

M 

0 

+ 

0 

+ 

M 

7. 

OH 

+ 

H 

H2 

+ 

0 

8. 

OH 

+ 

0 

-t- 

°2 

+ 

H 

9. 

OH 

+ 

H2 

->- 
-*- 

H20 

+ 

H 

10. 

OH 

+ 

OH 

-<- 

H20 

+ 

0 

11. 

CO 

+ 

OH 

co2 

+ 

H 

12. 

NO 

+ 

0 

<- 

°2 

+ 

N 

13. 

NO 

+ 

N 

-> 

N2 

+ 

0 

14. 

NO 

+ 

°2 

->- 
•t- 

N02 

+ 

0 
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TABLE    XI 


SELECTED    RATE    CONSTANTS 

FOR 

THE    ELEMENTARY    CHEMICAL    REACTIONS    FROM   NEWHALL     [10] 


Reaction 


Forward   Rate    Constant 


Backward   Rate 
Constant 


1 
2 
3 


5.4    x    1017   exp(-123,600/RT)  1.5    x   1016 


3.1xl015exp(-110,000/RT) 


7xl010     (150, 000)2 


,1/2 


(RT) 


12 


4.2x10  ,224,900 


J72 


5.4x10 


RT 
21 


exp(-74,000/RT) 


7.0    x    10 17    /  T 
14 


exp(-150,000/RT)  9.0    x    10 


-)exp(-224,900/RT)  6.1    x    10 


2.0    x    10 


14 


16 


6 

7 

8 

9 

10 

11 

12 

13 

14 


6xl013    11^000    exp(-ll,800/RT) 
12 


1.0    x    10 


12 


14 


1.4x10      exp (-6, 000/RT)  3 . 3x10      exp (-8 , 000/RT) 


13 


14 


5.5x10      exp (-1, 00  0/RT)  7. 2x10      exp (-16 , 90  0 /RT) 


13 


14 


6.2x10      exp (-6, 0  00/RT)  3.2x10      exp (-21, 100/RT) 


12 


13 


7.7x10      exp (-1, 0  0  0/RT)  8. 3x10      exp (-18 , 100/RT) 


12 


14 


7.1x10      exp (-7, 700/RT)  4.7x10      exp (-2 7,250 /RT) 

3.2xl09    exp(-39, 100/RT)       13 . 3xl09exp (-7 , 080/RT) 


1.55x10 


13 


7xl013exp(-75,500/RT) 


10ml/2 


1.8x10      Tx    exp (-47, 000/RT)  5.8x10 


10      ml/2 
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TABLE    XII 
RATE    CONSTANTS    AT    19  47°K 


Reaction 

Forward   Constant 

Backward 

Constant 

1 

7.1979    x    103 

1.5    x 

10" 

L6 

2 

1.3896    x    103 

3.5933 

X 

1014 

3 

3.4565    x    10"5 

9.0    x 

10 

L5 

4 

6.1003    x    10"10 

6.1   x 

10 

L4 

5 

1.36733    x    1010 

2.0    x 

10 

L6 

6 

1.037   x    102 

1.0    x 

10 

L4 

7 

2.9688    x    1011 

4.17295    : 

k    1011 

8 

4.247    x    1013 

9.1234 

X 

io12 

9 

1.3147    x    1013 

1.3693 

X 

io12 

10 

5.9461   x    1012 

7.7125 

X 

io11 

11 

9.702    x    1011 

4.1025 

X 

io11 

12 

2.542    x    108 

4.1536 

X 

io12 

13 

13 
1.55    x    10 

2.3418 

X 

IO5 

14 

4.205    x    106 

2.559 

X 

IO12 
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TABLE    XIII 

EQUILIBRIUM   MOLE    FRACTIONS   WITH   TOTAL   TEMPERATURE    =    2200 °K 
AND   TOTAL   PRESSURE    =    15    ATM    FROM   BANES,    ET    AL.  ,     [3] 


Throat 

Exit 

T 

1947°K 

1531°K 

P 

8.3    atm 

2.5  8    atm 

co2 

1.2919X10"1 

1.3073xl0-1 

H20 

1.2992xl0_1 

1.3048X10"1 

°2 

6.7092xl0"4 

5.3073xl0~5 

H2 

3.4679xl0"4 

3.6215xl0"5 

N2 

7.2896X10"1 

7.2981xl0_1 

N 

lxl0"12 

IxlO"12 

NO 

3.8768xl0~4 

2.4451xl0~5 

OH 

2.5755xl0~4 

1.2768xl0~5 

CO 

1.4098xl0~3 

1.0068xl0"4 

N02 

1.1142xl0"7 

2.8080xl0~9 

H 

7.9334xl0~6 

1.2395xl0"7 

0 

4.4208xl0~6 

3.8529xl0"6 

M 

28.876 

28.903 

P 

1.499xl0~3 

5.938xl0~4 
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TABLE    XIV 

EQUILIBRIUM   MOLE    FRACTIONS    WITH   TOTAL    TEMPERATURE    =    2200 °K 
AND   TOTAL    PRESSURE    =    5    ATM    FROM   BANES,    ET    AL.,     [3] 


Throat 

Exit 

T 

1947°K 

1531°K 

P 

2.7655    atm 

0.86    atm 

co2 

1.2845xl0_1 

1.3068xl0_1 

H20 

1.2963X10"1 

1.3045X10"1 

°2 

1.0156xl0"3 

7.6494xl0"5 

H2 

5.0574xl0"4 

5.0628xl0~5 

N2 

7.2857X10"1 

7.2979xl0_1 

N 

IxlO"12 

lxl0"12 

NO 

4.7613xl0"4 

2.9394xl0"5 

OH 

3.8269xl0"4 

1.8211xl0"5 

CO 

2.1936xl0~3 

1.7399xl0"4 

N02 

9.4140xl0"8 

2.4052xl0~9 

H 

1.7402xl0~5 

2.4370xl0"7 

0 

9.8841xl0"6 

7.6885xl0"8 

M 

28.863 

28.902 

P 

4.997xl0"4 

1.98xl0'4 
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TABLE    XV 
REACTIONS    RATES    OF   ELEMENTARY   EQUATIONS    AT    19  47°K 


Total   Pressure,    Atm 

15 

5 

Forward 

Backward 

Forward 

Backward 

1 

1.839xl0"6 

3.132xl0"6 

2.037xl0"7 

3.775xl0~7 

2 

9.48xl0"10 

2.3127x10 

1.534xl0~10 

4.114xl0"10 

3 

2.636xl0"17 

4.066xl0-16 

3.593xl0"18 

3.36xl0"16 

4 

8.748xl0"19 

1.02xl0~37 

-20 
9.703x10 

-24 

2.305x10 

5 

2.997xl0"6 

3.503xl0~6 

2.81xl0~7 

3.56xl0~7 

6 

1.369xl0~10 

1.997xl0~10 

2.30X10"11 

3.69xl0_11 

7 

1.637xl0~6 

1.727xl0"6 

5.92xl0~7 

6.25xl0"7 

8 

1.305xl0~4 

1.31xl0"4 

4.81xl0*5 

4.83xl0~5 

9 

3.169xl0~3 

3.809xl0~3 

7.62xl0"4 

9.26xl0"4 

10 

1.064xl0~3 

1.1955xl0"3 

2.61xl0~4 

2.28xl0~3 

11 

1.018xl0"3 

1.1348xl0"3 

2.44xl0"4 

2.75xl0"4 

12 

1.176xl0"9 

-12 
7.52x10 

3.58xl0~10 

-12 
1.26x10 

13 

1.620X10"11 

2.04xl0~9 

-12 
2.21x10 

5.05xl0"10 

14 

2.952xl0~9 

3.40xl0"9 

6.09xl0"10 

7.14xl0~10 
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TABLE    XVI 
EQUATIONS    FOR  TOTAL    REACTION    RATE    OF   EACH    SPECIE 


RC02  =  rllf"rllb 

^O  =  rlb"rlf+r9f"r9b+r10f"rllb 

R02  =  r6b"r6f+r8f~r8b+rl2f"r12b+r14b"r14f 

^2  =  r2b"r2f+r7f"r7b+r9b"r9f 

^2  =  r4b"r4f+r13f~r13b 

^0  r3b"r3f+r5f~r5b+rl2b"r12f+r13b~r13f+r14b"r14f 

ROH  rlf"rlb+r7b"r7f+r8b"r8f+r9b"r9f+2r10b"2r10f+rllb"rllf 

RCO  =  rllb        rllf 


^O 


«H 


5b      5f      14f      14b 
=      rlf-rlb+2r2f-2r2b+r7b-r7f+r8f-r8b+r9f-r9b+rllf-rllb 


R0  "      r3f    r3b+r5f"R5b+2r6f"2r6b+r7f~r7b+r8b'r8f+r10f"r10b 

+r12b"r12f+r13f"r13b+r14f~r14b 
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TABLE  XVII 

TOTAL  REACTION  RATE,  R,  AT  THE  THROAT 
WITH  STATIC  TEMPERATURE  =  19  47°K 


Species 


Total   Pressure 


15 


CO  -    1.168xl0-4  -  3.13xl0-5 

HO  -    7.7xl0~4  -  2.185xl0"3 

02  -    4.98xl0"7  -  1.695xl0~7 

H2  6.4xl0"4  1.63xl0~4 

N2  -   2.02xl0'9  -  5.03xl0~10 

NO  -    5.047xl0"7  -  7.445xl0~8 


OH  1.019xl0"3  4.24xl0"3 


CO  1.168xl0~4  3.13x10  5 


N02  5.056xl0"7  7.46x10    8 

H  -    7.59xl0~4  -    1.95xl0"4 


-    1.316xl0~4  -    2.022xl0"3 
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TABLE    XVIII 

CHARACTERISTIC    REACTION    TIMES    IN    THE    THROAT    AREA 
IN    ORDER   TO    REACH    EXIT    EQUILIBRIUM    CONDITIONS 
FOR   TOTAL    TEMPERATURE    =    2200 °K 

Total   Pressure 


15 


Species  Time 


C02  +   6.3x10    4 


H20  +   2.95xl0~5 


02  -   6.44xl0"~2 


H2  +    2.52xl0~5 

N2 

NO 


OH 


CO  +    6.26xl0"4 


N02  +    1.12xl0~5  +    3.62xl0_1 

H  -   5.34xl0~7  -   2.56xl0"2 


1.73xl0"6  -    1.40x10    3 


+ 

4. 

,32 

- 

3. 

,74x 

10 ' 

-2 

+ 

1. 

,25x 

10  " 

-5 

Transit  Time  5.7x10 


-5 


+ 

1.14 

- 

5.1xl0~3 

- 

9.59X101 

+ 

4.84xl0"2 

+ 

8.10xl03 

- 

1.04xl02 

+ 

1.49xl0"3 

+ 

2.00X101 
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TABLE  XIX 
FREEZING  CLASSIFICATION  OF  CHEMICAL  SPECIES 


Total 
Pressure 

Frozen 

-  ■  — 

Unknown 

i 

Equilibrium 

15 

°2 

H2° 

None 

N2 

co2 

NO 

H2 

H 

OH 

0 

CO 
N02 

5 

H20 

co2 

°2 
H2 

N2 

NO 

OH 

CO 

N02 

H 

0 

None 

None 
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APPENDIX      B 


Stator 


Exit    Guide    Vanes 


COMBUSTION 


CTfl 


CHAMBER 


y 


k 


Rotor 


^\r 


^-^-^ 


FIGURE    I        TURBINE    SCHEMATIC     AND     FLOW    STATION 
IDENTIFICATION     FOR    THE     VIBRATIONAL 

FREEZING    ANALYSIS 


59 


ENERGY 


Equilibrium 


Total 
Enthalpy7 


Rototiona 


Kinesis 


Energy 


ond  Translotional 
Energy 


p-Frozen   Fl 


low — m Relaxation  Zone 


3    AXIAL 
DISTANCE 


FIGURE  2 


ENERGY    DISTRIBUTION     IN     THE    TURBINE 
FOR   THE    VIBRATIONAL    FREEZING  ANALYSIS 
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Midspon  Streamline 


ENGINE 
AXIS        ~7 


/ 


FIGURE    3        STATOR    PASSAGE 
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LOCUS   OF 
MINIMA 

/ 
/ 


300  400  500  600  700 

CHANGE    IN     TOTAL   TEMPERATURE    ACROSS     THE 
TURBINE,      Ttz-    Tt'2,  °K 


FIGURE   4        PERCENTAGE      PERFORMANCE     LOST    DUE 
TO    VIBRATIONAL    FREEZING 
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15   ATM      N2-H20 


14  1.27  1.37 

AREA      RATIO,    A/A" 


.46 


FIGURE    5        LOCAL  VIBRATIONAL   RELAXATION    TIMES    FOR 
THE    STATOR    PASSAGE     IN    FIGURE       3 
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.96.. 


.92.. 


.90. 


.88 


.86 


th 


N2-N2    COLLISIONS 

MIDSPAN    STREAMLINE     LENGTH       6' 


.Equilibrium 


.14 


.27 


1.37 


1.46 


AREA    RATIO,    A/  A* 


.51 


FIGURE    6 


RATIO    OF   LOCAL    VIBRATIONAL    ENERGY     TO 
EQUILIBRIUM     THROAT    VIBRATIONAL     ENERGY 
AS  A  FUNCTION     OF   AREA    RATIO    FOR 
PURE    NITROGEN 
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96-- 


.94-- 


*    i 


vth 

.921 


.90.. 


.88.. 


.86 


SYMBOL 


MIDSPAN  STREAMLINE    COLLISION 
LENGTH  PARTNERS 

6"  N2-N2 

2"  N2~N2 

2"  N2  -N2 

N2-H20 


-EQUILIBRIUM 


.14  1.27  1.37 

AREA    RATIO,     A/ A* 


.46 


1.51 


FIGURE  7  INFLUENCE  OF  STATOR  SIZE  AND  COLLISION 
PARTNER  ON  THE  RATIO  OF  LOCAL  VIBRATIONAL 
ENERGY  TO  EQUILIBRIUM  THROAT  VIBRATIONAL 
ENERGY  AS   A   FUNCTION   OF  AREA    RATIO 
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COMPLETE       EQUILIBRIUM 


ENTHALPY 


1 1 

y/*\    - 

pt 

2 

3 

3, 

^-^2 

ENTROPY 


CHEMICALLY 

FROZEN 

I, 

i 

ENTHALPY 

J^ 

- 

^A 

ENTROPY 


FIGURE     8 


EXPANSION     PROCESSES      USED    IN 
ANALYSIS    OF    CHEMICAL   FREEZING 
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700 
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-T, 


K 


FIGURE    9 


PERCENTAGE    PERFORMANCE      LOST     AS      A 
RESULT    OF     CHEMICAL      FREEZING 
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APPENDIX  C 

An  approximate  solution  to  equation  (21)  may  be  obtained 
using  the  finite  difference  method.  Equation  (21)  is  put  in 
the  form  as  in  Hurle,  et  al.,  [9]. 

[1  -  (Ax/2Xn)] 

E (Tv} n+1   =   [1  +  (Ax/2X  ~)  ]    E (Tv} n 

n+1 

A    Ax/2^  E(T)n+l    E(T)n 


[1  +  (Ax/2X  TtTj   aTt:      X 
'       n+1      n+1       n 

where  n  and  n+1  refer  to  two  adjacent  locations  along  the 
nozzle  axis  separated  by  distance  Ax,  and  X  and  T  are  evalu- 
ated for  equilibrium  flow.   X  is  the  characteristic  relaxation 
length  and  X  =  Ut ,    where  U  is  the  local  gas  velocity  and  t  is 
relaxation  time. 
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APPENDIX  D 

This  is  the  calculation  of  the  conversion  constant  for 
use  in  the  rate  equations  so  that  mole  fractions  may  be 
utilized  directly. 

The  units  of  the  rate  constants  in  Table  XI  vary  with 
the  type  of  elementary  equation;  bi-molecular  or  tri- 
molecular.   The  normal  units  for  concentrations  are  moles 
-  cm   and  the  units  of  X  are  moles  of  specie  i  per  mole  of 
mixture.   From  the  ideal  gas  law 

P  V  =  n  R  T 
Multiplying  by  moles  of  specie  i  over  moles  of  specie  i 

P  V  X.   =   n.  R  T 
i       1 


or 


n.      P  X. 


V       R  T 

-3 
in  units  of  moles  -  cm   ;  the  desired  units. 

Using  the  gas  constant  equal  to 

r   =   82.06  Cn\  a^ 
mole  °K 

a  conversion  factor  for  use  at  a  specified  temperature  and 

mixture  pressure  can  be  obtained.   The  mole  fractions  are  each 

multiplied  by  this  factor  and  are  then  ready  for  use  in  the 

specie  rate  equations. 
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